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The HI and ZEUS Collaborations have performed new next-to-leading order QCD analyses to deter- 
mine the parton density functions of the proton. QCD fits are performed using inclusive neutral and 
charged current deep inelastic scattering cross sections from HERA-I. The fits include a full treat- 
ment of experimental systematic uncertainies, taking into account point-to-point correlations. The 
extracted parton densities are in agreement with those from global fits. Since HERA inclusive data 
provide no direct information on the high- a; gluon, an independent fit has been performed by the 
ZEUS Collaboration in which the inclusive DIS cross sections are supplemented by jet data from deep 
inelastic scattering and photoproduction. The determination of the gluon distribution is significantly 
improved in such fits, allowing a competitive extraction of the strong coupling, a B (Mi). 



1 Introduction 

The kinematics of lepton-proton scattering 
are described in terms of the Bjorken scal- 
ing variable, x, the negative invariant mass 
squared of the exchanged vector boson, Q 2 , 
and the fraction of energy transferred from 
the lepton to the hadronic system, y. 

At Leading Order (LO), in the elec- 
troweak interaction, the double differential 
cross section for Neutral Current (NC) Deep 
Inelastic Scattering (DIS) can be written in 
terms of structure functions, 



d 2 a± c 



[Y+F 2 - y 2 F L =F y_xF 3 ] 



dxdQ 2 

where Y± = 1 ± (1 — y) 2 . Similarly, for the 
Charged Current (CC) process, 



[F + F 2 CC± - y 2 F™± T Y.xF 3 



cc±i 



dxdQ 2 

where the CC structure functions depend on 
the charge of the incoming lepton. The struc- 
ture functions are directly related to the par- 
ton density functions (PDFs) and the Q 2 de- 
pendence, or scaling violation, is predicted in 
perturbative QCD. 

Conventionally, QCD analyses use the 
formalism of the next-to- leading order (NLO) 
DGLAP evolution equations^. These equa- 
tions yield the PDFs at all values of Q 2 , pro- 
vided they are input as functions of x at some 



starting scale Qq. The x parameterisation at 
Qq is usually chosen to be of the general form, 

xf(x) = Ax b (l-x) c P{x) 

where P(x) is a polynominal function. While 
A controls the normalisation, the parameters 
b and c are sensitive to the low and high- a; 
regions, respectively. 



2 HERA-Only QCD Analysis 

The PDFs of the proton are usually deter- 
mined in global fits to both colliding beam 
data from HERA, as well as to DIS data 
from fixed target machines. In such analy- 
ses, the high precision NC cross sections from 
HERA are crucial in determining the low- a; 
sea and gluon distributions, while the fixed 
target data provide most of the information 
on the high- a; sea and gluon, as well as on 
the valence quark distributions. In the global 
fits, the most important inputs for the deter- 
mination of the valence PDFs have been the 
vFe and the fiD fixed target data. However, 
these can suffer from large uncertainties due 
to heavy target correction^ 

In the present analyses, the HI and 
ZEUS Collaborations have performed QCD 
fits to only HERA data, thus elminating 
any uncertainties due to heavy target cor- 
rections. The full set of inclusive DIS data 



from HERA-I has been used, which covers 
a large range in (x,Q 2 ). Since NC and CC 
DIS data provide no information on the high- 
er gluon, an independent fit has been per- 
formed by ZEUS in which the inclusive data 
are supplemented by jet cross sections in 
DIS (Q 2 » 1 GeV 2 ) and photoproduction 
(Q 2 ~ GeV 2 ) . 

The PDFs from the HI and ZEUS analy- 
ses are presented with full accounting for un- 
certainties from correlated systematic uncer- 
tainties, as well as from statistical and uncor- 
rected sources. 

2.1 The HI PDF 2000 Analysis 

A full description of the HI PDF 2000 anal- 
ysis is given elsewherJSl Here, only a sum- 
mary is presented. The inclusive DIS dat£0 
used in the analysis, span the kinematic range 
8 x 10~ 5 < x < 0.65 and 1.5 < Q 2 < 
30000 GeV 2 . To ensure the applicability 
of perturbative QCD, an additional cut of 
Q 2 > 3.5 GeV 2 is imposed on the data in- 
cluded in the fit. The parameterisation of the 
PDFs, at the starting scale of Qo = 4 GeV 2 , 
is of the general form, 

xf{x) = Ax b {l - x) c (l + ex + fx 2 + gx 3 ). 

The PDFs that are parameterised are the 
U = u + c (total up-type), D = d + s (total 
down-type), U, D and g and the number of 
terms in the polynomial is chosen separately 
for each PDF by searching for \ 2 saturation. 
The fit has a total of 10 free parameters and 
the analysis is performed in the Zero Mass 
Scheme, which is appropriate for high Q 2 . 

The resulting parton densities, at a scale 
of Q 2 = 4 GeV 2 , are shown in Fig.Q] The u- 
and d-valence distributions are constructed 
from u v = U—U and d v — D—D respectively. 
The dark shaded band shows the total exper- 
imental uncertainty, while the light shaded 
band shows the model uncertainty, which in- 
cludes contributions from a variation of the 
input scale Qq, the minimum Q 2 cut on the 
data, the charm and strange fractions, the 
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Figure 1. The parton densities from HI PDF 2000 at 
Q 2 = 4 GeV 2 . The dark shaded band is the experi- 
mental uncertainty and the light shaded band is the 
model uncertainty. The solid line shows the results 
of the Hl+BCDMS fit. 

quark masses and the value of a s . The lat- 
ter results in the largest contribution to the 
uncertainty on the gluon. The fit provides 
a tight constraint on the total up-type and 
down-type quark densities at low-x. 

A cross-check to the HI PDF 2000 fit has 
also been performed, in which the HI inclu- 
sive NC and CC data are supplemented with 
precise fixed target data from BCDM$3 The 
resulting PDFs are shown by the solid line 
in Fig. □ The Hl+BCDMS fit is consis- 
tent with the results of HI PDF 2000. The 
HI PDFs are also found to be compatible 
with those from the global fits of MRSTPand 
CTECP (see Fig.©. 

2.2 The ZEUS-Only PDF Fit 

Full details of the ZEUS analysis are given 
elsewhere^. Here, the main features are 
briefly described. The inclusive DIS dat£0 
used in the fit, span the kinematic range 

6.3 x 10~ 5 < x < 0.65 and 2.7 < Q 2 < 



30000 GeV 2 . To reduce higher twist effects, 
an additional cut of W 2 > 20 GeV 2 is im- 
posed. The form of the parameterisation is, 

xf(x) = Ax b (l-x) c {l + ex) 

at Ql — 7 GeV 2 . No advantage in the x 2 
is found to result from using more complex 
polynomial forms. The PDFs that are pa- 
rameterised are the u v (u- valence), d v (d- 
valence), S (total sea), g and A = (d — u) a . 

Since the data from HERAT are less 
precise than fixed target data in the high- a; 
regime, the high-x sea and gluon are not well 
constrained in fits using only HERA data. In 
the ZEUS-Only fit, the high- a; sea and gluon 
parameters, cs and c g , are fixed to the val- 
ues extracted fro m th e previously published 
ZEUS-S global frJES. The constraints im- 
posed lead to a total of 10 free parameters. 
The fit is performed in the Roberts-Thorne^l 
variable flavour number scheme. 

The ZEUS-Only PDFs are summarised 
in Fig. |3 and are compared to those of the 
ZEUS-S global fit. The central values of 
MRST and CTEQ are also indicated. The 
high- a; valence is not quite as well constrained 
in the ZEUS-Only fit as in the global fits (c.f. 
ZEUS-S). However, they are becoming com- 
petitive, particularly for the less well known 
(i-valence distribution. Furthermore, these 
PDFs are free from the potentially large un- 
certainties from heavy target corrections. 

At low-x, the sea and gluon PDFs are 
just as well determined as in the correspond- 
ing distributions of the global fits, since 
HERA data are crucial in determining these 
distributions for all fits. At high-x, the 
ZEUS-Only PDFs have the same uncertain- 
ties as the global fits because of the con- 
straints that have been applied to the high- a; 
sea and gluon. 

"Since there is no information on the shape of the 
d—u distribution in fits to HERA data alone, this dis- 
tribution has its shape fixed consistent with the Drcll- 
Yan datsll^I and its normalisation cons isten t with the 
size of the Gottfried sum-rule violation^ 
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Figure 2. The PDF distributions from the ZEUS- 
Only fit compared to those of the ZEUS-S global fit, 
with full uncertainties. The central values of MRST 
and CTEQ are also shown. 

The experimental uncertainties represent 
the most significant source of uncertainty 
on these distributions. Variation of analysis 
choice, such as the value of Qq, the minimum 
Q 2 of data entering the fit and changing the 
form of the parameterisation at Qg, do not 
produce a large model uncertainty. 

2.3 The ZEUS- Jets PDF 

In QCD fits incorporating only inclusive NC 
and CC DIS data, the high- a; gluon is con- 
strained by the momentum sum rule only. 

The ZEUS Collaboration have performed 
a new fit, in which the inclusive DIS data are 
supplemented by cross sections from inclu sive 
jet DIjJMI and two-jet photoproductior^l 
Both sets of jet data are directly sensitive to 
the gluon in the range 0.01 < x < 0.1. 

A full NLO calculation for jet cross sec- 
tions is very time-consuming. Therefore, 
NLO programs 6 have been used only initially 

b For the inclusive DIS jets, the program DISENtED 
was used, while for the two-jet photoproduction, the 
program of Frixione and Ridolfi - ' was used. 



For Publisher's use 




Figure 3. The PDF distributions from the ZEUS-Jets 
and HI PDF 2000 fits, with full uncertainties. The 
central values of MRST and CTEQ are also shown. 

to produce grids of weights, giving the pertur- 
batively calculable part of the cross section. 
The predictions for the jet cross sections can 
then be reconstructed using, 



E 



da; a™{n R ) /„(£, fiF)ca,n(pF, to) 



where f a is the PDF for parton a as a func- 
tion of momentum fraction £ and scale fiF, 
and the c a , n are the perturbatively calculable 
kernels. The cross sections calculated accord- 
ing to this equation reproduce the real NLO 
predictions to better than 1%. 

Recall that for the ZEUS-Only fit 
(Sec. 12.211 . the high- a; gluon parameter, c g , is 
constrained to be consistent with the ZEUS- 
S global fit. However, the extra information 
from the jet data allows this parameter to be 
freed, giving 11 free parameters in total. 

The PDF distributions extracted from 
the ZEUS-Jets fit are compared in Fig. [21 to 
those of HI PDF 2000. The central values of 
MRST and CTEQ are also shown. All fits 
are compatible within uncertainties. 

The jet data are expected to have most 
impact on the gluon distribution. Figure 0] 



shows the gluon PDFs for fits without c (left) 
and with (right) jet data. The inclusion of 
jet data provides a significant improvement 
to the constraint on the gluon at medium-to- 
high- a;, which persists to high scales. 

Extraction of a s (M|) 

The value of a s (M§) has been determined 
from the ZEUS-JETS fit by treating it as 
an additional free parameter. The value ex- 
tracted is, 

a s (M§) = 0.1183 ±0.0027 (experimental) 
± 0.0008 (model) ± 0.0050 (scale) 

where the experimental uncertainty arises 
from both uncorrelated and correlated 
sources and the model uncertainty includes 
contributions from varying the value of Qq, 
changing the form of the input parameteri- 
sation and varying the cuts on the data in- 
cluded in the fit. The uncertainty in a s (M^), 
which usually comes from the correlation to 
the PDF shapes, is automatically included in 
the experimental uncertainties. The largest 
uncertainty comes from varying the scale /in, 
suggesting that NNLO QCD analyses could 
provide very precise extractions of a s in the 
future. 

3 Summary 

The HI and ZEUS Collaborations have per- 
formed NLO QCD DGLAP analyses of the 
HERA-I data to extract the parton density 
functions of the proton. Since HERA is a 
proton-only target machine, these analyses 
avoid any uncertainties due to heavy tar- 
get corrections, which have led to large sys- 
tematic uncertainties in global fits. In addi- 
tion, since the data used in the fits are from 
only one experiment, uncertainties which can 
arise from combining systcmatics from diffcr- 



c This fit is labelled "ZEUS-O" since it includes the 
same data as the ZEUS-Only fit described in Sec. 12. 21 
However the parameter constraints are the same as 
for the ZEUS-JETS fit i.e. free c a . 
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Figure 4. The gluon distribution for fits without 
(left) and with (right) jet data included. The cross- 
hatched band show the uncorrelated uncertainty 
while the shaded band shows the total uncertainty. 

ent experiments are also reduced. Jet data 
have been added in an independent fit, giv- 
ing a significant improvement in the knowl- 
edge of the gluon PDF at medium-to-high- x 
which persists to high scales. The extra con- 
straint on the gluon has led to a competitve 
extraction of a s (M§) from only HERA data. 

Acknowledgments 

I would like to thank my colleagues in the 
ZEUS and HI QCD Fitting Groups. This 
work was supported by PPARC. 

References 

1. V.N. Gribov and L.N. Lipatov, Sov. J. 
Nucl. Phys 20, 94 (1975); G. Altarelli 
and G. Parisi, Nucl. Phys. B126, 298 
(1977); Yu.L. Dokshitzer, Sov. Phys. 
JETP 46, 641 (1977) 

2. A.M. Cooper-Sarkar et al, J. Phys G25, 
1387 (1999) 

3. HI Coll., C. Adloff et al, Eur. Phys. 
J. C13, 609 (2000); HI Coll., C. Adloff 
et al, Eur. Phys. J. C19, 269 (2001); 
HI Coll., C. Adloff et al., Eur. Phys. J. 



C21, 33 (2001); HI Coll., C. Adloff et 
al, Eur. Phys. J. C30, 1 (2003) 

4. HI Coll., Eur. Phys. J. C30 (2003) 

5. BCDMS Coll., A.C. Benvenuti et al, 
Phys. Lett. B223, 485 (1989) 

6. A.D. Martin et al Eur. Phys. J. C, 73 
(2002); A.D. Martin et al Eur. Phys. 
J. C, 455 (2002) 

7. J. Pumplin et al, JHEP 0207 , 012 
(2002) 

8. ZEUS Coll., S. Chekanov et al, Eur. 
Phys. J. C21, 443 (2001); ZEUS Coll., 
J. Breitwig et al, Eur. Phys. J. C12, 
411 (2000); ZEUS Coll., S. Chekanov et 
al, Eur. Phys. J. C28, 175 (2003); 
ZEUS Coll., S. Chekanov et al, Phys. 
Lett. B539, 197 (2002); ZEUS Coll., S. 
Chekanov et al, Preprint DESY-03-214 
| |hep-ex/0401003| ), DESY, 2004, Subm. 
to Phys. Rev. D.; ZEUS Coll., S. 
Chekanov et al, Eur. Phys. J. C32, 
16 (2003) 

9. ZEUS Coll., "An NLO QCD analysis 
of inclusive cross section data and jet 
production data from the ZEUS exper- 
iment at HERA-I", submitted to the 
32nd International Conference on High 
Energy Physics, 16-22 August 2004, Bei- 
jing, China 

10. R.S. Towell et al, Phys. Rev. D64, 
052002 (2002) 

11. P. Amaudruz et al, Phys. Lett. B295, 
159 (1992); P. Amaudruz et al, Phys. 
Rev. Lett. 66, 2712 (1991) 

12. S. Chekanov et al, Phys. Rev. D67, 
012007 (2003) 

13. R.G. Roberts and R.S. Thornc, Phys. 
Rev. D57, 6871 (1998) 

14. ZEUS Coll. S. Chekanov et al, Phys. 
Lett. B547, 164 (2002) 

15. ZEUS Coll. S. Chekanov et al, Eur. 
Phys. J. C23, 615 (2002) 

16. S. Frixione and M.H. Seymour, Nucl. 
Phys. B485, 291 (1997) 

17. S. Frixione and G. Ridolfi, Nucl. Phys. 
B507, 315 (1997) 



